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Abstract 

We show that the phase difference of two overlapping ground state Bose- 
Einstein condensates can effect the optical spontaneous emission rate of ex- 
cited atoms. Depending on the phase difference the atom stimulated sponta- 
neous emission rate can vary between zero and the rate corresponding to all 
the ground state atoms in a single condensate. Besides giving control over 
spontaneous emission this provides an optical method for detecting the con- 
densate phase difference. It differs from previous methods in that no light 
fields are applied. Instead the light is spontaneously emitted when excited 

atoms make a transition into either condensate. 
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Spontaneous symmetry breaking is a fundamental physical phenomenon. In particular, 
breaking of the global phase symmetry explains much of the interesting physics of Bose- 
Einstein condensates [^. We show that this spontaneously broken symmetry can effect the 
rate of spontaneous emission from excited atoms. Hence we are able to link two fundamental 
physical phenomena: spontaneous symmetry breaking and spontaneous emission. 

Other physics which effects the spontaneous emission rate are; atom stimulation by a 
single condensate |^ , non-free space electromagnetic boundary conditions [§] , due to a mirror 
for example, and superradiant type dipole correlations [^. The latter occurs when the atoms 
are much closer than the wavelength of the emitted light. Our effect is different to these and 
happens when the excited atoms can decay into either of two final states which are Bose 
condensed. Its physical origin is the indistinguishability of the final states, which produces 
phase dependent interference terms in the transition probability. 

Two groups have reported producing two dilute gas Bose-Einstein condensates in the 
same trap. Spatially overlapping condensates of the \F = l,m = —1) and |F = 2,m = 
2) hyperfine spin states of rubidium-87 have been produced by sympathetic cooling [Q. 
Separated condensates of the same state of sodium-23 have been produced using a far blue 
detuned laser sheet to divide a magnetic trap |]5[. Our system requires two overlapping 
condensates into which an excited state can make a spontaneous transition. 

Bose stimulation of transitions by atoms has been studied for both electronic |@] and 
nuclear transitions. It was shown to be experimentally feasible for electronic transitions 
0. We apply the same idea to a system with two ground state condensates instead of one, 
and find that the atom stimulated emission rate depends on the condensate phase difference. 

This phase difference has been observed directly by destructive optical imaging of the 
interference fringes in the density of overlapping condensates 0. Previously proposed meth- 
ods for optically detecting this phase difference require driving Raman transitions between 
the ground states. The excited state was created only through coherent excitation of the 
ground states and hence was not independent of the ground states. Javanainen has shown 
that prompt amplification of one Raman beam by light-stimulated scattering is a signature 
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of a condensate phase difference 0. This is because the condensate phase difference be- 
haves hke a Raman optical coherence. Imamoglu and Kennedy have studied hght scattering 
from two spatially separated condensates. They showed that the optical coherence allows 
spontaneous Raman scattering to generate coherent light, as well as allowing the complete 
elimination of Raman scattering 0. These two schemes were treated qualitatively, in con- 
trast to the fully quantitative treatment we shall give. Ruostekoski and Walls found that the 
relative peak heights of the characteristic two-peaked incoherent light scattering spectrum 



depend on the relative condensate phase Our method is related to these three insofar as 
it results from the macroscopic quantum coherence between the two condensates. However 
no external light is applied and no Raman transitions occur. The light is emitted by the 
decay of independently produced excited atoms. A schematic experimental setup is shown 
in Fig. 1. 

We consider transitions from an excited state labeled by e to ground states labeled by 
gi and g2. The excited state is separately prepared and then launched towards the ground 
state condensates. The ground states may differ in their internal quantum numbers and/or 
in their external wavefunctions. Both situations have been achieved experimentally [^Q. 
We neglect radiation reaction and multiple scattering effects |rT|,|T2|. This requires that 
the system be optically thin at the frequency of the spontaneously emitted light We 



show that the recoil shift of the emitted photon makes this possible, provided the natural 
linewidth is sufficiently small. Hence we have in mind the decay of a metastable excited 
state. 

We consider spontaneous emission from excited atoms. "Spontaneous" specifies that 
the emission is not stimulated by light, although it may be stimulated by atoms. The 
spontaneous emission is from the excited state atom field ijje to ground state atom fields i^g^ . 
The positive frequency part of the electric field operator at position r is given in terms of 

as 



the atom fields by Javanainen and Ruostekoski [13 



E+(r) = E+(r) + ^ / rfV K(d,; r, r')^J(r')4(r'), (1) 
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where E^(r) is the free field, which we assume to be the vacuum so that (E]^(r) •E^(r)) = 0. 
Rotation at the transition frequency Q has been removed from the excited atom and electric 
fields, so that they are slowly varying [|l^. We ignore excited state depletion and ground 
state growth. The integration kernel K(d; r, r') generates the field at r due to a dipole 



moment d at r' [|TJ]. We shall only be concerned with electromagnetic waves in the far field 
so that the dominant part of the kernel is that proportional to |r — r'|^^, 

/N k"^ expfi/clrl) , ., , 

K(dg; r, r') ^ ^^r^ exp{-ikr' ■ n)X>„ (2) 

'±tt6q |r| 

where dg is the dipole moment of the transition g ^ e, and /c = |k| where k is the wavevector 
of the emitted light. The far field assumption, in the form |r| ^ |r'|, has been used to replace 
|r — r'l in the denominator by |r| and to expand the exponential exp(iA;|r — r'|) ^ exp(iA;|r| — 
ikr'-n), where n = r/|r|. We have also approximated the propagation direction (r— r')/|r— r'| 
by n, and hence introduced the components of the dipole moments perpendicular to the 
propagation direction n, 

X>3 = (nxdg)xn = dg-n(n-dg). (3) 

With these approximations the expectation value of the field intensity, with respect to 
the atom field states and a vacuum free field, is 



J(r)=2c£o(E-(r)-E+(r)) 

= OiE/ ^'r'/ dh" exp[-ik{T" - r') • n] 



X 
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X»*-X»,,(^t(r')^^(r')^t,(r")4(r")), (4) 



where k = ck^ / {^Tt'^SQ). As usual we consider the atom field operators to be expanded in 
complete sets of orthonormal mode functions {0j(r)} which are adapted to the system under 
consideration 'iP{y) = J2i d,i(f)i{T), where the a.j are the usual mode annihilation operators. The 
ground state modes into which condensation occurs are typically solutions to the Gross- 



Pitaevskii equation |15 
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We assume a sufficiently low temperature that stimulation by the non-condensed atoms 
can be ignored. In a JILA experiment up to 80% of the trapped atoms have been measured 



to be condensed, at a temperature of about 140 nK ||I6[. Since the ratio of non-condensed 
to condensed atoms is then 1/4 this implies that spontaneous emission stimulated by non- 
condensed atoms can be reduced, at least, to 1/4 of that stimulated by the condensates. 
This would be a constant background emission. It would add to the unstimulated free space 
spontaneous emission rate 7, which is always present. Furthermore, in their two condensate 
experiment the JILA group reports cooling proceeding until "we can no longer see any 
noncondensed atoms" @]. 

We assume the Bose condensed modes may be represented by coherent states 
\\/N exp(i^)). The coherent state phase 9 embodies the broken global phase symmetry. 
Since Hamiltonians always contain pairs of atom creation and annihilation field operators 
the phase is unmeasureable. However phase differences between condensates are measure- 
able, although experimentally they are not yet controllable and hence they vary from shot to 
shot. The coherent states have a mean atom number of N and an atom number uncertainty 
of ^/N. We further assume that all condensates are initially uncorrelated. The assumption 
that the excited state is condensed is not essential, however it does simplify our treatment. 
The atom field operator expectation value determining the intensity Eq.(^ is then 



(^t(r')^,(r')^J,(r")4(r")) = N^^/n^' exp(z5e,,,)0:(r')0.(r')0;'(r")0e(r"), (5) 

where 69gg/ = Og — Ogi is the ground state condensate phase difference. The atom stimulated 
intensity can then be written as 

= i^^e \Y.Ng\Cg'Dg? + 2^Ng,Ng,Re {exp(z5^,,,jX>; ■ , (6) 

where Re{} means the real part, and we have defined Franck-Condon factors, or the overlaps 
of the excited and ground modes taking into account the photon recoil, 

Cg = j dh'exp{-tkr' ■ n)0;(r')0e(r'). (7) 
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The last term in the expression for the intensity Eq.(|^) depends on the condensate phase 
difference 59g^g^ provided the Franck- Condon factors are non-zero and T>*^ ■ "Dg^ ^ 0. That 
is provided there is some field polarization that both dipole moments can generate at r. 
Since the electric field is a vector field this is necessary for interference to be possible. This 
dependence of the spontaneously emitted intensity on the condensate phase difference could 
in principle be used to measure it. 

We now give a physical explanation based on the quantum mechanical transition ampli- 
tudes. First we recall that if the final states of two transitions are in principle distinguishable 
then the individual transition probabilities must be added to get the final state probabil- 
ity. However if the final states are indistinguishable then the transition amplitudes must 
be added, before taking the squared modulus to get the final state probability. The squar- 
ing produces interference terms in the final state probability. The intensity radiated in the 
propagation direction n is proportional to the transition probability for producing a pho- 
ton propagating in that direction. The first terms in Eq.(|^) correspond to the individual 
probabilities for the transitions into each ground state. The last term corresponds to the 
product of amplitudes for each transition and hence represents quantum mechanical inter- 
ference between them. Interference occurs because the condition T>*^ ■ "Dg^ ^ ensures that 
the particular transition which produces photons with certain polarizations cannot be dis- 
tinguished. A general treatment of the electric field polarization dependence on the atomic 
level scheme is given by Javanainen and Ruostekoski ||13|| . 

We next estimate the mode overlaps Cg. We use gaussians for the ground mode functions 

0, = {2'nl')-'"eM-{^TDe,f/{Al% (8) 

where is the one dimensional variance of the probability density 0^, is the unit vector 
in the z direction, and we have allowed for possible displacement by the distance ±D in 
this direction. For the excited atoms we use the plane wave mode (pe = V^^^''^ exp(ifcex). 
This represents atoms propagating in the positive x direction with momentum hke- V is 
the mode volume. The plane wave usefully models modes which are much larger than the 
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ground modes, so for purely calculational reasons we require (2/)^/V^ ^ 1. We assume that 
the excited state atoms are completely independent of the ground state atoms. They would 
be prepared in a separate part of the experiment. Evaluating the integral Eq.(|^) we find 

Cg = (87r)3/4 / _ J exp[-Z2(fc2 + - 2kke cost?)] x 

exp [±ikD sin t? cos </)] , (9) 

where t? is the angle between the light emission direction n and the excited state propagation 
direction e^,., and (f is the polar angle about e^.. The first exponential results from momentum 
conservation while the second is a phase modulation due to displacement of the condensate 
from the coordinate origin. 

Since the mode functions {(pi} are a complete ortho normal set the fraction of emissions 
into the solid angle dQ with final atomic mode (pg is jC^P dQ. Choosing the photon and 
excited state momenta equal k = ke, which gives the optimal overlap, and using the realistic 
values l^/V = 0.01, and Ik = 100 gives 

|CJ2 ~ l.3exp[-4 X 10^(1 -cost?)] 

^ 1.3exp[-2 X 10^^% (t? < 1). (10) 

This mode overlap is largest for emission into a cone about the excited state propagation 
direction p|. For the atom stimulated emission rate to exceed the non-stimulated rate by a 
factor of one hundred, at an emission angle of t? = 0.02 radians, requires about 100/|Cgp ~ 
2.3 X 10^ atoms in the stimulating condensate. At t? = 0.03 radians the stimulated emission 
is negligible. 

A condition for the validity of our analysis is that the spontaneously emitted light is 
sufficiently far off resonance to avoid cooperative effects and multiple scattering |T^. As- 
suming k = fee, the emitted photon is detuned from resonance by the recoil frequency 
ujR = hk'^ / (2m), where m is the atomic mass. This detuning must be at least of the order of 
the condensate collective linewidth, which has been estimated to be SNg'y /{2{lk)'^), where 7 
is the transition's free space natural linewidth [0,|l^. The condition is then 



O^R > - 



(11) 
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With Ng = 10^ and Ik = 100 this becomes ujr > I5O7. A typical recoil frequency of 
ujji = 10^/s requires the natural linewith 7 < 600/s. We note that the total transition rate 
is proportional to N^. which is not constrained. 

We next consider some particular examples. Let the excited state have magnetic quantum 
number m. It can decay by emission of (a") or (cr+) circularly polarised photons, or by 
emission of a linearly polarised photon. The respective final states having magnetic quantum 
numbers m + 1, m — 1 and m. We assume that transitions into the m ground state, which is 
not condensed, can be ignored. This is justified when transitions into the condensed mil 
modes are dominant because of Bose stimulation. 

The total radiated power is found by integrating Eq.(|D, for the intensity radiated in 
direction n, over all directions. Let z be the quantization axis and let the excited atoms 
propagate in the positive x direction, see Fig. 1. We use the circular dipole moments 
d± = ^d±(e,j. ± iey)/^/2, where d± = \d±\ are the magnitudes of the dipole moments, and 
the e^/y are unit vectors in the x/y directions. In the z direction these dipoles radiate 
photons with opposite circular polarization. In the x/y plane, however, they radiate indis- 
tinguishable linearly polarized photons, which makes interference possible, see Fig. 1. Using 
the previously calculated mode overlaps Eq.(|^) the total atom stimulated radiated power is 
found to be 



The integral is over the spherical surface of radius |r| centered on the origin. We have defined 
the coefficients 




(12) 




B 



cosh a 



+ 



(a^ — 1) sinh a 



(13) 
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where a = For the case of complete overlap D = the remaining coefficient has the 

simple form 

^ ^sinha 

Cd=o = 3B- 4 . (14) 

a 

In the limit of ground state mode functions much larger than the wavelength of light, a ^ 1, 
the power reduces to 

Pn=o = N,N,K{87r)'/'-^^d\l - cos69+^), (15) 

where for simplicity we have assumed here, and hereafter, that = = Ng, = d_ = d, 
and k = kg. This expression shows that the total spontaneously emitted power may vary 
substantially with the condensate phase difference 56'+_. We note that the excited state 
decay rate is proportional to the total emitted power. This means that the depletion of the 
excited state, as measured by direct detection of excited state atoms, could also be used to 
detect the phase difference. 

The power Eg . (|T5|) can be scaled to an emission rate by dividing by the photon energy 
hck. Using the free space spontaneous emission rate 7 = d'^k^/{3iTeoh) and the numerical 
values preceding Eq.([TO|) gives the total atom stimulated emission rate 

7..». = 7iV.iV,-L_L_(^-j (l-cos^^,_) 

= -fNeNg{2A X 10~^)(1 - cos50+_). (16) 



For cos 69^ = this is one hundred times the free space emission rate •yN^ when A^, 

4.2 X 10^. 

We next consider spatially separated condensates, that is with _D 7^ in Eq.(|). The 
coefficient C of the interference term in Eq.(|l2]) for the total radiated power is then 

C - V^t^ (^f + , (17) 

where is the modified Bessel function of the ffist kind of order n. The terms in the series 
proportional to result from the overlap of the gaussian wavef unctions. These terms are 
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dominant for wavefunction widths much bigger than the wavelength, a ^ 1. A graph of 
this coefficient as a function of the ratio of condensate width to separation is shown in Fig. 
2, for kl = Att. This shows that the effect of the interference on the total radiated power 
disappears once the condensates are separated. 

The physical situation is related to the Young's two slit interference experiment with 
slits of width 21 and separation 2D. In the one dimensional case the interference pattern 
is the product of a sine function due to the single slit diffraction and a cosine due to the 
interference 

/(m) = |^^^^^cos(A:/m + ^,)|', (18) 

where m = sin ■(9, ■(9 is the angle of observation of the interference pattern, and 6'^ is the phase 
difference between the light waves at each slit. This phase difference determines the relative 
positions of the interference pattern and the diffraction pattern. For the two slit experiment 
D > I. However overlapping condensates have D < I. Consider the effect of varying the 
source phase difference 9d in this case. Since the interference maxima are further apart than 
the sine function zeros, Eq.(|TB|) shows that the interference maxima can shift in and out 
of under the maxima of the sine function as 9d varies. This is the physical origin of the 
modulation of the total power by the condensate phase difference. When D > I this does 
not happen because the interference maxima are closer together. Hence when one maximum 
shifts out another shifts in to replace it and the total power is approximately constant. The 
overall effect is well approximated by considering the effect on the total power to arise from 
the overlapping portions of the condensates. 

As a final example we consider the polarization dependence of the radiation for another 
geometry. The positive frequency part of the electric field with polarization Cj, perpendicular 
to n, is the scalar product E+(r) ■ e*, with E+(r) given by Eq.dTj). The intensity of the 
polarised field is given by similar expressions to Eq.(^ and Eq.(^ with T^g replaced by 
'Dg ■ e* = dg ■ e* . Let the quantization, observation, and excited atom propagation directions 
all be z. With this geometry the transitions are distinguishable by the polarization of the 
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emitted photon and the total intensity does not depend on the relative condensate phase. 
However the intensity in particular linear polarizations does depend on the relative phase, 
as we now show. Define the linear polarization vectors e/3 = e^; cos (3 + ey sin (3. Then the 
intensity in polarization is 

1/3 = NeNg^^ + - Re{C:C_ exp[z(2/3 + 5^+-)]}) . (19) 

This intensity is polarization dependent due to the final term, which also depends on the 
ground state condensate phase difference. This polarization dependence is diagnostic of 
broken symmetry in the form of a relative phase between the ground state condensates. 
When the relative phases of the transition dipoles of different atoms are independent there 
is no such polarization dependence. Independence is impossible for condensed atoms, which 
are all in the same state. However for non-condensed atoms independence arises through 
collisions or other dephasing perturbations. 

We summarise our results by listing measurable quantities that may depend on the 
ground state condensate phase difference: the intensities of particular electromagnetic field 
modes, the total radiated power, and the total transition rate. The physical origin of these 
effects is the indistinguishability of the final states. 

An experiment relying on atom stimulated emission into two condensates will be more 
difficult than the corresponding single condensate experiment. However our numerical esti- 
mates, and the recent experimental interest in two condensates in the same trap give 
cause for optimism concerning the feasibility of observing the effect of spontaneously broken 
phase symmetry on spontaneous emission. 

This work was supported by the Marsden Fund of the Royal Society of New Zealand, 
The University of Auckland Research Fund and The New Zealand Lottery Grants Board. 
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FIGURES 

FIG. 1. Schematic diagram of an experiment investigating the condensate phase dependence 
of spontaneously emitted light. The solid ellipses represent the ground state condensates and the 
dashed ellipse represents the excited mode. The excited mode has momentum hk in the x direction. 
The dipole moments are quantized in the z direction and are shown as arrows and directed circles. 
Note that when viewed from the x direction both dipoles appear to oscillate linearly. 

FIG. 2. Graph of the interference coefficient C, Eg. (|T7|) , as a function of condensate separation 
D. Parameters are kl = Att and a = 4(47r)^, corresponding to condensates larger than a wavelength 
in size. 
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